Out-of-equilibrium phase transitions in the HMF model: a closer look 
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We provide a detailed discussion of out-of-equilibrium phase transitions in the Hamiltonian Mean 
Field (HMF) model in the framework of Lynden-Bell's statistical theory of the Vlasov equation. 
For two-levels initial conditions, the caloric curve only depends on the initial value /o of the 

distribution function. We evidence different regions in the parameter space where the nature of phase 
transitions between magnetized and non-magnetized states changes: (i) for /o > 0.10965, the system 
displays a second order phase transition; (ii) for 0.109497 < fo < 0.10965, the system displays a 
second order phase transition and a first order phase transition; (iii) for 0.10947 < fa < 0.109497, the 
system displays two second order phase transitions; (iv) for fo < 0.10947, there is no phase transition. 
The passage from a first order to a second order phase transition corresponds to a tricritical point. 
The sudden appearance of two second order phase transitions from nothing corresponds to a second 
order azeotropy. This is associated with a phenomenon of phase reentrance. When metastable states 
are taken into account, the problem becomes even richer. In particular, we find a new situation 
of phase reentrance. We consider both microcanonical and canonical ensembles and report the 
existence of a tiny region of ensembles inequivalence. We also explain why the use of the initial 
magnetization Mq as an external parameter, instead of the phase level fo, may lead to inconsistencies 
in the thermodynamical analysis. 

PACS numbers: 



I. INTRODUCTION 



Systems with long-range interaetions have recently 
been the object of an intense activity [iHl]- These 
systems are numerous in nature and concern different 
discipUnes such as astrophysics (galaxies) 0-01: two- 
dimensional turbulence (vortices) [8|-[lC|. biology (chemo- 
taxis) [ll| , plasma physics [T^ - [l4 | and modern technolo- 
gies such as Free Electron Lasers (FEL) p^| - [T7| . In ad- 
dition, their study is interesting at a conceptual level 
because it obUges to go back to the foundations of sta- 
tistical mechanics and kinetic theory P. [T8l - l20| . Indeed, 
systems with long-range interactions exhibit a number 
of unusual features that are not present in systems with 
short-range interactions. For example, their equilibrium 
statistical mechanics is marked by the existence of spa- 
tially inhomogeneous equilibrium states [l| , unusual ther- 
modynamic limits |2l| - [23j , inequivalence of statistical en- 
sembles 0, 0, [111 ; negative specific heats [2^ , vari- 
ous kinds of phase transitions 0, [l^l etc. Their dynam- 
ics is also very interesting because these systems can be 
found in long-lived quasi stationary states (QSS) that 
arc different from Boltzmann equilibrium states. These 
QSSs can be interpreted as stable steady states of the 
Vlasov equation which governs the evolution of the sys- 
tem for sufficie ntly "short" times before correlations have 
developed [H, HI]. In fact, for systems with long-range 
interactions, the coUisional relaxation time towards the 
Boltzmann distribution increases rapidly with the num- 
ber of particles N and diverges at the thermodynamic 
limit N +00 0, [H, HO, Hi] . Therefore, the domain 
of validity of the Vlasov equation is huge and the QSSs 



have very long lifetimes. In many cases, they are the 
only observable structures in a long-range system, so that 
they arc often more physically relevant than the Boltz- 
mann equilibrium state itself. A question that naturally 
emerges is whether one can predict the QSS actually 
reached by the system. This is not an easy task since 
the Vlasov equation admits an infinite number of stable 
steady states in which the system can be trapped [11] . In 
a seminal paper, Lynden-BcU (29j proposed to determine 
the QSS eventually reached by the system by develop- 
ing a statistical mechanics of the Vlasov equation. To 
that purpose, he introduced the notions of phase mixing, 
violent relaxation and coarse-grained distributions. He 
obtained the most probable distribution by maximizing 
a Boltzmann-type entropy while conserving all the con- 
straints of the Vlasov equation (in particular the infinite 
class of Casimirs). By definition, this "most mixed state" 
is the statistical equilibrium state of the Vlasov equation 
(at a coarse-grained scale). Whether or not the system 
truly reaches this equilibrium state relies on an assump- 
tion of ergodicity and efficient mixing. This ergodicity 
assumption is not always fulfilled in the process of vio- 
lent relaxation and the Lynden-Bell prediction may fail. 
In that case, the QSS can be another stable steady state 
of the Vlasov equation that is incompletely mixed. This 
is referred to as incomplete relaxation (see, e.g. [soj . for 
discussion and further references). In case of incomplete 
relaxation, the prediction of the QSS is very difficult, 
and presumably impossible. Nevertheless, in many cases, 
the Lynden-Bell approach gives a fine first order predic- 
tion of the achieved QSS and allows one to predict out- 
of-cquilibrium phase transitions between different types 
of structures that can be compared with direct Simula- 
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tions or experiments. Before addressing this problem in 
a specific situation, namely the Hamiltonian Mean Field 
(HMF) model [HI, [s^l, let us first briefly review the suc- 
cesses and the weaknesses of the Lynden-Bell approach. 

Lynden-BcU's statistical theory of violent relaxation 
was elaborated in the context of 3D stellar systems. Un- 
fortunately, this is the worse situation for its practical 
application. Indeed, the predicted distribution function 
has infinite mass (the spatial density decreases at large 
distances like In other words, this means that there 

is no entropy maximum for a stellar system in an infinite 
domain [g^. This is a clear evidence of the fact that 
galaxies have necessarily reached a state of incomplete 
violent relaxation. In fact, the Lynden-Bell theory is able 
to explain the isothermal core of elliptical galaxies with- 
out recourse to collisions that operate on a much longer 
timescale (of the order of the Chandrasekhar relaxation 
time (33!). This is usually recognized as a major suc- 
cess of the theory. Unfortunately, it fails at predicting 
the structure of the halo whose velocity distribution is 
anisotropic and whose spatial density decreases like 
j28j . Models of incomplete violent relaxation have been 
elaborated by Bertin & Stiavelli [s^l, Stiavelli & Bertin 
(ssj and Hjorth & Madsen [3^. These models are able 
to reproduce the dc Vaucouleurs law of elliptical galax- 
ies and provide a very good agreement with numerical 
simulations up to nine orders of magnitude [37j . An- 
other possibility to describe incomplete relaxation is to 
develop a kinetic theory of violent relaxation in order 
to understand what limits mixing [38l - l40| . The idea is 
that, in case of incomplete relaxation (non-ergodicity), 
the prediction of the QSS is impossible without consid- 
ering the dynamics (soj . Finally, in more academic stud- 



4l| . one can confine the system within an artificial 



spherical box and assume a complete relaxation inside 
the box. Since the Lynden-Bell distribution is similar to 
the Fermi-Dirac statistics (in the two-levels approxima- 
tion), the problem is mathematically equivalent to the 
study of a gas of self-gravitating fermions in a box. This 
theoretical problem has been studied in detail by Chava- 
nis [l^j. The caloric curve P{E) displays a rich variety of 
microcanonical and canonical phase transitions (zcroth 
and first order) between gaseous (non degenerate) and 
condensed (degenerate) states, depending on the value 
of a degeneracy parameter related to the initial distribu- 
tion function /o in the Lynden-Bell theory. In particular, 
there exists two critical points in the phase diagram, one 
in each ensemble, at which the phase transitions are sup- 
pressed. For details about these phase transitions, and 
for an extended bibliography, we refer to the review 0. 
The Lynden-Bell prediction has also been tested in ID 
and 2D gravity [43l |43 | where the infinite mass problem 
does not arise j45| . However, it is found again that relax- 
ation is incomplete and that the Lynden-Bell prediction 
fails nil- Finally, Arad & Lynden-Bell have shown 
that the theory itself presents some inconsistencies aris- 
ing from its non-transitive nature. These negative results 
have led many astrophysicists to the conclusion that the 



Lynden-Bell theory does not work in practice [28|. 

A similar statistical theory has been developed by 
Miller [l^], and independently by Robert & Sommeria 
[48j . in 2D turbulence in order to explain the robust- 
ness of long-lived vortices in astrophysical and geophys- 
ical fiows (a notorious example being Jupiter's great red 
spot). Large-scale vortices arc interpreted as quasi sta- 
tionary states of the 2D Euler equation in the same way 
that galaxies are quasi stationary states of the Vlasov 
equation (see [l^, IS] for a discussion of the numerous 
analogies between the statistical mechanics and the ki- 
netic theory of 2D vortices and stellar systems). Miller- 
Robert-Sommeria (MRS) developed a statistical theory 
of the 2D Euler equation in order to predict the most 
probable state achieved by the system. Although situa- 
tions of incomplete relaxation have also been evidenced 
in 2D turbulence [sol - is^ . the MRS theory has met a lot 
of success. For example, it is able to account for geom- 
etry induced phase transitions between monopoles and 

toles as we change the aspect ratio of the domain [53l - 
. Phase transitions and bifurcations between different 
types of fiows have also been studied in [58l - l60| . On the 
other hand, when applied to geophysical and astrophys- 
ical flows, the MRS theory is able to account for the 
structure and the organization of large-scale fiows such 
as jovian jets and vortices f6ll-l64| and Fofonoff flows in 
oceanic basins [HI, |65| . This theory has also been applied 
to more complicated situations such as 2D magnctohy- 
drodynamics (MHD) [66l . [67| and axisymmctric flows (the 
celebrated von Karman flow) [g^ ■ 

A toy model of systems with long-range interactions, 
called the Hamiltonian Mean Field (HMF) model, has 
been introduced in statistical physics [lH, 112 and exten- 
sively studied |3] . It can be viewed as a XY spin system 
with inflnite range interactions or as a one dimensional 
model of particles moving on a ring and interacting via a 
long-range potential truncated to one Fourier mode (co- 
sine potential). In that second interpretation, it shares 
many analogies with self-gravitating systems (sU . [3^ |69| 
but is much simpler to study since it avoids difficulties 
linked with the singular nature of the gravitational po- 
tential at the origin and the absence of a natural confine- 
ment [l-l^ • The observation of quasi stationary states in 
the HMF model [13, [z3| was a surprise in the community 
of statistical mechanics working on systems with long- 
range interactions. It was recognized early that these 
QSSs are out-of-equilibrium structures and that they are 
non-Boltzmannian. They were first interprete d [TOl in 
terms of Tsallis generalized thermodynamics |7l| with 
the argument that the system is noncxtensive so that 
Boltzmann statistical mechanics is not applicable. Later, 
inspired by analogies with stellar systems and 2D vor- 
tices reported in [lOj] , different groups started to interpret 
these QSSs in terms of stable steady states of the Vlasov 
equation and statistical eq uilibrium states in the sense 
of Lynden-Bell [H, IH 1^. Chavanis (tJ studied out- 
of-equilibrium phase transitions in the HMF model by 
analogy with similar studies in astrophysics and hydro- 
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dynamics j4lU53ll and obtained a phase diagram in the 
(/o, E) plane [901 between magnetized (M 7^ 0) and non- 
magnetized (M = 0) states. These regions are separated 
by a critical line Ed fo) that marks the domain of stabil- 
ity of the homogeneous phase. This critical line displays 
a turning point at ((/o)*,-B*) ~ (0.10947,0.608) leading 
to a phenomenon of phase reentrance (as we reduce the 
energy, the homogeneous phase is successively stable, un- 
stable and stable again). Antoniazzi et al. |75| studied 



the validity of the Lynden-Bell prediction by performing 
careful comparisons with direct A^-body simulations at 
E = 0.69 and found a good agreement for initial magne- 
tizations Mo < {Mo)critiE) ~ 0.897 lea ding to spatially 
homo gen eous Lynden-Bell distributions [9 If . Antoniazzi 
et al. [78| obtained a phase diagram in the (Mq, E) plane 
and showed that the system exhibits first and second or- 
der phase transitions separated by a tricritical point. Fi- 
nally, Antoniazzi et al. [79| performed numerical simu- 
lations of the Vlasov equation and found a good agree- 
ment with direct A^-body simulations and Lynden-Bell's 
prediction for the explored range of parameters. A syn- 
thesis of these results was published in [s^]. In this pa- 
per, a more detailed discussion of phase transitions in 
the {fo,E) plane was given, showing the lines of first 
and second order phase transitions and the domains of 
metastability. On the other hand, a comparison between 
the phase diagrams in the (/o, E) and (Mq, E) planes was 
made. It was stated, without rigorous justification, that 
the tricritical point in the {Mo,E) plane corresponds to 
the turning point of the critical line Edfo) in the (/o, E) 
plane, i.e. the point where the phase reentrance starts. 
These results were confronted to numerical simulations 
by Staniscia et al. (sTj . These simulations confirmed the 
existence of a reentrant phase in the very narrow region 
predicted by the theory [tII but also showed discrepan- 
cies with the Lynden-Bell prediction (such as an addi- 
tional reentrant phase and a persistence of magnetized 
states in the a priori non-magnetized region) that were 
interpreted as a result of incomplete relaxation. Stanis- 
cia et al. [8l| also determined the physical caloric curve 
Pkin{E), where /3fci„ = l/Tkm is the inverse kinetic tem- 
perature, in the region of the phase diagram displaying 
first and second order phase transitions, and reported 
the existence of a region of negative kinetic specific heat 
Ckin = dE/dTkin < 0. In a recent paper [s^ . the thermo- 
dynamical caloric curve I3{E) was determined in the same 
region of parameters and it was shown that the thermo- 
dynamical specific heat C = dE/dT is always positive, 
even in the region where the kinetic specific heat is neg- 
ative. In particular, it is argued that the ensembles are 
equivalent although the experimentally measured specific 
heat is negative |92| . 

These various results show that the description of out- 
of-equilibrium phase transitions in the HMF model is 
extremely rich and subtle. In this paper, we describe 
in more detail the phase transitions between magnetized 
and non-magnetized states in the {fo,E) plane. In par- 
ticular, we plot the series of equilibria f3{E) for differ- 



ent values of fo and determine the caloric curve cor- 
responding to fully stable states. This completes and 
illustrates our previous study HH where only the final 
phase diagram was reported. We evidence different re- 
gions in the parameter space where the nature of phase 
transitions changes: (i) for fo > {fo)t — 0.10965, the 
system displays a second order phase transition; (ii) for 
(/o)i ~ 0.109497 < fo < (/o)t - 0.10965, the system 
displays a second order phase transition and a first or- 
der phase transition; (iii) for (/o)* — 0.10947 < fo < 
(/o)i ~ 0.109497, the system displays two second order 
phase transitions; (iv) for fo < (fo)* — 0.10947, there is 
no phase transition. The passage from a first order phase 
transition to a second order phase transition corresponds 
to a tricritical point. The sudden appearance of two sec- 
ond order phase transitions from nothing corresponds to 
a second order azcotropy. This is associated with a phe- 
nomenon of phase reentrance. When we take into account 
metastable states, the description is even richer and seven 
regions must be considered (see Sec. IIIip . In particular, 
we find a new situation of phase reentrance. We also 
stress two unexpected results that were not reported (or 
incorrectly reported) in previous works: (i) Contrary to 
what is stated in [S^l, there exists a region of ensembles 
inequivalence but it concerns an extremel y n arrow range 
of parameters so that the conclusions of |82l | are not al- 
tered; (ii) the tricritical point separating second and first 
order phase transitions does not exactly coincide with 
the turning point of the stability line Ec{fo), contrary to 
what is stated in [s^l, but is slightly different. Again, 
the difference is small so that the main results of previ- 
ous works are not affected. However, this slight difference 
leads to an even richer variety of phase transitions. We 
may be fascinated by the fact that so many things hap- 
pen in such a very narrow range of parameters (typically 
(/o),„ ~ 0.1075 < fo < {fo)c ^ 0.11253954) although 
fo can take a priori any value between and +00! Fi- 
nally, we make clear in this paper (see Sec. |ll| that the 
relevant control parameters associated with the Lynden- 
Bell theory are {fo, E) [t^I while the use of the variables 
{Mo,E) [zl, [7^ may lead to physical inconsistencies in 
the thermodynamical analysis. 



II. THE LYNDEN-BELL THEORY AND THE 
CHOICE OF THE CONTROL PARAMETERS 

The HMF model [3l|, H^l , which shares many similar- 
ities with gravitational and charged sheet models, de- 
scribes the one-dimensional motion of A^ particles of unit 
mass moving on a unit circle and coupled through a mean 
field cosine interaction. The system Hamiltonian reads 

N N 

1=1 i,j=l 

where 6i represents the angle that particle i makes with 
an axis of reference and Vi stands for its velocity. The 
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1/N factor in front of the potential energy corresponds 
to the Kac prescription to make the system extensive 
and justify the vahdity of the mean field approximation 
in the limit N — >■ +00. The relevant order parameter 
is the magnetization defined as M = mij/N where 
m.; = (cos6'i, sin^i). In the N — > +00 limit, the time 
evolution of the one body distribution function f(9,v,t) 
is governed by the Vlasov equation 



^ + .^-(A4[/]sm^ 



My[f] cos 9)^=0, (2) 



where M^f] = J f {9, v,t) cosO dOdv and My[f] = 
J f(9, V, t) sin 9 d9dv arc the two components of the mag- 
netization. 

The statistical theory of the Vlasov equation, intro- 
duced by Lynden-Bel l [29|. h as been reviewed in several 
papers ||, [ll [zl, III EoTili so that we shall here only 
recall the main lines that are important to understand 
the sequel. We assume that the initial distribution func- 
tion takes only to values f{9,v,t = 0) G {0,/o}. For 
example, it can be made of one or several patches of 
uniform distribution f{6,v,0) = fo surrounded by "vac- 
uum" f{9,v,Q) = 0. We note that the number and the 
shape of these patches can be completely arbitrary. For 
such initial conditions, the quantities conserved by the 
Vlasov equation are: (i) the value fo of the initial dis- 
tribution; (ii) the normalization A4 = J f d9dv = 1; (iii) 
the energy E = ^ J fv'^ d9dv + 5(1 - Af^). The fine- 
grained distribution function f{9,v,t) is stirred in phase 
space but conserves its two values fo and at any time, 
i.e. f{9,v,t) £ {0, /o} Vt. However, as time goes on, 
the two levels values fo and become more and more in- 
termingled as a result of a mixing process (filamentation) 
in phase space. The coarse-grained distribution f{9, w, t), 
which can be viewed as a local average of the fine-grained 
distribution function, takes values intermediate between 
and fo, i.e. < f{9, v, t) < fo. It is expected to achieve 
a steady state f{9, v) as a result of violent relaxation on a 
relatively short timescale (a few dynamical times). This 
corresponds to the QSS observed in the simulations. The 
most probable, or most mixed state, is obtained by max- 
imizing the Lynden-Bell entropy 
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d9dv. 



(3) 



while conserving E and ^4 (for a given value of /q). 
This determines the statistical equilibrium state of the 
Vlasov equation. Note that the whole theory relies on 
an assumption of ergodicity, i.e. efficient mixing. Our 
aim here is not to determine the range of validity of the 
Lynden-Bell theory, so that we shall assume that this as- 
sumption is fulfilled (see, e.g. [tJ, for a discussion of 
incomplete relaxation in the HMF model). We are led 
therefore to considering the maximization problem 



m^x{S[f]\E[f] = E,M[f] = l}. 
f 



(4) 



for a given value of fo- The critical points of cancel- 
ing the first order variations of the constrained entropy, 
are given by the variational principle 



6S - PSE - aSM = 0, 



(5) 



where and a are Lagrange multipliers. This yields the 
Lynden-Bell distribution 
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where e{9,v) = v"^ /2 - M^[f^B\cos9 - My[/iB]sin0 
is the individual energy. In the two-levels approxima- 
tion, the Lynden-Bell distribution is formally identical 
to the Fermi-Dirac statistics Note that T = = 
{dS/dE)~^ is the thermodynamical temperature. Since 
the distribution function ^ is non-Boltzmannian, the 
thermodynamical temperature differs from the classical 
kinetic temperature Tkin = C fv'^ d9dv. This point has 
been studied specifically in [S^ . 

The maximization problem (j4]) corresponds to the mi- 
crocanonical ensemble (MCE). Since the Lynden-Bell 
theory is based on the Vlasov equation that describes 
an isolated system, the microcanonical ensemble is the 
relevant ensemble to consider (the energy is fixed). We 
can, however, formally define a canonical ensemble. We 
introduce the free energy functional J[f] = S[f] — /3E[f] 
[93^ and consider the maximization problem 



max{j[/]|A^[/] = l} 
/ 



(7) 



for a given value of fo. The maximization problems (|4]) 
and ([7]) have the same critical points since the variational 
principle 



5J - aSM = 



(8) 



returns Eq. ([5|) (recall that /? is fixed in the canonical 
ensemble). In addition, it can be shown at a general 
level [2^ that a solution of the canonical problem ([7]) is 
always a solution of the more constrained dual micro- 
canonical problem but that the reciprocal is wrong 
in case of ensembles inequivalence [oj]. Therefore, even 
if the canonical ensemble is not physically justified in 
the context of Lyndcn-BcU's theory of violent relaxation, 
it provides nevertheless a sufficient condition of micro- 
canonical thermodynamical stability. It is therefore use- 
ful in that respect. In addition, it is interesting on a 
conceptual point of view to study possible inequivalence 
between microcanonical and canonical ensembles. There- 
fore, we shall study in this paper the two maximization 
problems ([4]) and ([7]), while emphasis and illustrations 
will be given for the more physical microcanonical case. 

Before that, let us recall general notions that will be 
useful in the sequel (for an extended account, see e.g. 
Q). For a given value of fo, the series of equilibria is 
the curve P{E) containing all the critical points of ^ 
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or ([7]) (as we have seen, they are the same). The sta- 
ble part of this curve, in each ensemble, gives the corre- 
sponding caloric curve. In MCE, the control parameter 
is the energy and the stable states are maxima of en- 
tropy S at fixed energy and normalization. This defines 
the microcanonical caloric curve P{E). In CE, the con- 
trol parameter is the inverse temperature and the stable 
states are maxima of free energy J at fixed normaliza- 
tion. This defines the canonical caloric curve E{/3). The 
strict caloric curve contains only fully stable states (S) 
that are global entropy maxima at fixed energy and nor- 
malization in MCE or global free energy maxima at fixed 
normalization in CE. The physical caloric curve contains 
fully stable and metastable states (M), that are local en- 
tropy maxima at fixed energy and normalization in MCE 
or local free energy maxima at fixed normalization in 
CE. The unstable states (U), that are minima or sad- 
dle points of the thcrmodynamical potential, must be 
rejected. Note that for systems with long-range inter- 
actions, metastable states can have very long lifetimes 
so that they are very important in practice. By study- 
ing the caloric curve P{E) for a given value of /q, we 
can describe phase transitions. Microcanonical first or- 
der phase transition are marked by the discontinuity of 
the inverse temperature /?(£') at some energy Ef. This 
corresponds to a discontinuity of the first derivative of 
entropy S'{E) = I3{E) at Et in the energy vs entropy 
curve. There can exist metastable branches around Et 
that possibly end at microcanonical spinodal points. Mi- 
crocanonical second order phase transitions are marked 
by the discontinuity of P'{E) at some energy Ec. This 
corresponds to a discontinuity of the second derivative of 
entropy S"{E) ~ P'{E) at Ec. Similarly, canonical first 
order phase transitions are marked by the discontinuity 
of energy E{f3) at some inverse temperature f3f This 
corresponds to a discontinuity of the first derivative of 
free energy J'{E) = —E{/3) at /3t in the inverse temper- 
ature vs free energy curve. There can exist metastable 
branches around Pt that possibly end at canonical spin- 
odal points. Canonical second order phase transitions 
are marked by the discontinuity of E' {j3) at some inverse 
temperature j5c. This corresponds to a discontinuity of 
the second derivatives of free energy J"(/3) = —E'{(3) at 
Pc- Finally, by varying the external parameter /o, we can 
describe changes from different kinds of phase transitions 
at some critical values of /o and plot the corresponding 
phase diagrams (/o, E) and (/o, (3) in microcanonical and 
canonical ensembles. This is the programm that we shall 
follow in this paper. 

We emphasize that these general results are valid for 
the caloric curve /?(£') where /3 is the inverse thcrmo- 
dynamical temperature, not the inverse kinetic temper- 
ature. In particular, the thcrmodynamical specific heat 
C = dE / dT is always positive in the canonical ensemble 
while the kinetic specific heat Ckin = dE/dTkm can be 
positive or negative in the canonical ensemble. This has 
been illustrated in (s^ . 

A last comment is in order. If we consider a watcrbag 



initial condition in which f{9,v,t = 0) = /o in the rect- 
angle [0,nin,Omax] X [-Vmin, V„iaT] and f{0,V,t = 0) = 

outside, it seems convenient to take as control param- 
eters the initial magnetization AIq and the energy E as 
done in [tI,!?^. Indeed, the specification of these param- 
eters determines /o — (piE, Mq) and E and thus allows 
to compute the corresponding Lyndcn-BcU state. There- 
fore, it seems that the choice of the control parameters 
{E,Mo) or (Ej/q) is just a question of commodity. In 
fact, this is not the case, and we would like to point out 
some difficulties in taking {E, Mq) as control parameters 
in the thcrmodynamical analysis: 

(i) The control parameters {E,Mo) are less general 
than {E, fo) because they assume that the initial con- 
dition is a waterbag distribution, whereas the control 
parameters {E, /□) are valid for any initial distribution 
with two levels, whatever the number of patches and their 
shape. They allow therefore to describe a wider class of 
situations. 

(ii) The variables {E, Mq) may lead to redundancies 
because there may exist two (or more) couples {E, Mg^^) 

(2) 

and {E,Mq ) that correspond to the same {E,fQ) and, 
consequently, to the same Lyndcn-BcU state (recall that 
the Lyndcn-Bell prediction only depends on E and fo) 
[H. This has been illustrated in [80,llli. 

(iii) More importantly, the use of Mq as an external pa- 
rameter (instead of fo) leads to physical inconsistencies 
in the thcrmodynamical analysis. Indeed, if we work in 
terms of the variables {E, Mq), the initial value of the dis- 
tribution fo becomes a function fo = 4'{E, Mo) of these 
variables. As a result, the Lyndcn-BcU entropy functional 
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dedv, (9) 



depends not only of the external parameter Mq but also 
on the energy E. This is clearly a very unconventional 
situation. Indeed, if we want to apply the standard re- 
sults recalled above, the entropic functional can depend 
on an external parameter but it cannot explicitly depend 
on the energy. Therefore, these general results [2J] are 
not valid for functionals of the form In particular, 
the "improper" caloric curve I3{E) at fixed Mq can dis- 
play a region of negative specific heat while the proper 
caloric curve I3{E) at fixed fo does not. This is exempli- 
fied in Fig. 1(b) of [79j where the entropy versus energy 
is plotted at fixed Mo- This curve has a convex dip (re- 
vealing a negative specific heat region), while the curve 
S(E) at fixed fn has no convex dip and the ensembles 
are equivalent |82| . 

Finally, in the other contexts where the Lynden-Bell 
theory has been applied [4l|, \E^, [58l - [60l [63j , the control 
parameters that have been taken are E and /q. It is 
therefore important to describe the phase transitions in 
terms of these parameters as initiated in [t^ . 
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0.105 0.11 0.115 0.12 0.125 



FIG. 1: Phase diagram in the microcanonical ensemble. We 
have indicated several representative points. The minimum 
accessible energy of the homogeneous phase, correspond- 
ing to a waterbag (or Fermi) distribution, is Emin{fo) = 
l/(967r^/o^) + 1/2 [Zi, [131. We have also indicated the mini- 
mum accessible energy (for homogeneous and inhomogeneous 
phases) in the case where the initial condition is a rectangu- 
lar waterbag initial condition. It is equal to Emin{fo) when 
^< 0.12135... and to £m/jv(,/o) when fo > 0.12135... (see 
[8ll | for details). The curves Et{fo) and Em{fo) have been 
continued "by hand" (due to numerical problems) and may 
not be correct for small energies (see Appendix [X)) . 



III. DESCRIPTION OF CALORIC CURVES 
AND PHASE TRANSITIONS 

A. Phase diagrams 

In Fig. [1] we reproduce the microcanonical phase di- 
agram obtained in [8i[. In Fig. [21 we enlarge this dia- 
gram close to the turning point of energy ((/o)*,£'*) — 
(0.10947,0.608) in order to show that its structure is 
more complicated than previously thought. Similarly, in 
Figs. |3] and IH we plot the canonical phase diagram and 
its enlargement close to the turning point of temperature 
((/o)*,/?*) ^ (0.10947, 118). These phase diagrams show 
that we must consider different regions where the nature 
of phase transitions changes. 

In the following sections, we plot the series of equilib- 
ria 13 {E) in seven characteristic regions of the phase di- 
agram and describe the corresponding phase transitions. 
The branches (S) correspond to fully stable states, the 
branches (M) correspond to metastable states and the 
branches (U) correspond to unstable states. At the end of 
each subsection, we summarize the nature of phase tran- 
sitions in the corresponding region by considering only 
fully stable states. 

We will find that the microcanonical and canonical 
phase transitions are very similar. In fact, the ensembles 
differ only in a very small range of parameters. There- 
fore, we will essentially focus on the microcanonical en- 
semble and only mention the canonical ensemble in case 
of ensembles inequivalence. 
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FIG. 2: Enlargement of the phase diagram in the microcanon- 
ical ensemble. 
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Phase diagram in 


the canonical ensemble. 



B. Region 5 

In Figs. [S][7]we plot the series of equilibria in Region 5 
corresponding to fo > (/o)e where (/o)c = l/(27r-\/2) ~ 
0.11253954 (see Fig. [J) QHH- Specifically, we consider 
fo = 0.1130. 

The homogeneous phase exists at any accessible en- 
ergy. It is fully stable for E > and unstable for 
E < Ec- The inhomogeneous phase exists for E < Ec- It 
has a higher entropy (see Fig. |6|) than the homogeneous 
phase and it is fully stable. Therefore, the microcanoni- 
cal caloric curve displays a second order phase transition 
between homogeneous and inhomogeneous states marked 
by the discontinuity of f3'{E) a.t E ~ Ec- In the entropic 
curve of Fig. [HI this corresponds to a discontinuity of the 
second derivative S"{E) = j3'{E) at E = E^. The mag- 
netization passes from M = for E > Ec to M ^ fox 
E < Ec but remains continuous at the transition (see Fig. 
[7]). The discussion is similar in the canonical ensemble. 

Region 5: (i) in the MCE, there exists a second order 
phase transition at Ec- (ii) In the CE, there exists a 
second order phase transition at /3c. The ensembles are 
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FIG. 4: Enlargement of the phase diagram in the canonical 
ensemble. 



FIG. 6: Entropy versus energy in Region 5. 
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FIG. 5: Series of equilibria in Region 5 displaying microcanon- 
ical and canonical second order phase transitions at Ec and 
/3c respectively. 



equivalent. 
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In Figs. [SlfTOl we plot the series of equilibria in Region 
4 corresponding to (/o)t < /o < (/o)c where (/o)t ~ 
0.10965 and (/o)c 0.11253954 (see Fig. [1]). Specifically, 
we consider /o = 0.1110. 

The homogeneous phase exists at any accessible en- 
ergy. It is fully stable for E > unstable for 
E^^ < E < e'-P and metastablc for E < . A first 
inhomogcncous phase exists for E < Ec^ '' . It has a higher 
entropy than the homogeneous phase and it is fully sta- 
ble (see Fig. [9]). Therefore, the microcanonical caloric 
curve displays a second order phase transition between 
homogeneous and inhomogeneous states marked by the 
discontinuity of S"{E) = 13' {E) at E = E^ . The mag- 
netization passes from M = for E > i^c^' to M 7^ crgy. 



1.58 0.59 0.6 0.61 0.62 0.63 0.64 0.65 



FIG. 7: Magnetization versus energy in Region 5. 



(see Fig. [TU)) . A second inhomogeneous phase exists for 

(2) (2) 

E < Ec ■ It appears precisely at the energy Ec at 
which the homogeneous phase becomes metastable. It 
has a lower entropy S than the homogeneous phase and 
the first inhomogeneous phase (sec Fig. |9]) and it is un- 
stable. This branch is clearly visible on the magnetiza- 
tion curve (see Fig. [T0| . The discussion is similar in the 
canonical ensemble. 

Region 4: (i) in the MCE, there exists a second order 
phase transition at Ec^\ (ii) In the CE, there exists a 



second order phase transition at (3, 
equivalent. 



(1) 



The ensembles are 



for E < E^^^ but remains continuous at the transition 



D. Region 3-c 

In Figs. [TTIfTBl we plot the series of equilibria in Region 
3-c corresponding to (/o)2 < /o < (/o)t where (/o)2 ^ 
0.109519 and (/o)t =i 0.10965 (see Fig. ^. Specifically, 
we consider /o = 0.10963. 

The homogeneous phase exists at any accessible en- 
unstable for 
E < e''^ and 
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It is fully stable for E > E,_ 
metastablc for Et < 



eP <e<e[' 
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FIG. 8: Series of equilibria in Region 4 displaying microcanon- 
ical and canonical second order phase transitions at Ei^"^ and 
/3c^' respectively. At E^^^ and I3i^\ a second inhomogeneous 
phase appears (but is unstable) while the homogeneous phase 
becomes metastable. 



FIG. 11: Series of equilibria in Region 3-c. It displays mi- 
crocanonical and canonical second order phase transitions at 
and jSl^^ respectively. It also displays microcanonical 
and canonical first order phase transitions at Et and /3t re- 
spectively. 
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FIG. 9: Entropy versus energy in Region 4. 
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FIG. 10: Magnetization versus energy in Region 4. The un- 
stable inhomogeneous phase is clearly visible. 



fully stable for E < Et- A first inhomogeneous phase ex- 
ists for E < E^J-K It is fully stable for Et < E < E^'^^ and 
metastable for E < Ef. Indeed, it has a higher entropy 
than the homogeneous phase for Et < E < E^^^ and a 
lower entropy for E < Et- Therefore, the microcanoni- 
cal caloric curve displays a second order phase transition 
between homogeneous and inhomogeneous states marked 
by the discontinuity of /3'{E) = S"{E) at E = E^'^^ (see 
Fig. [T^ and a first order phase transition between ho- 
mogeneous and inhomogeneous states marked by the dis- 
continuity of (3{E) = S'{E) at E = Et (see Fig. [H]). 
The magnetization of the fully stable branch passes from 
M = to A/ 7^ at = i?!^' but remains continuous, 
and it passes from M ^ to M = at E = Et by being 
discontinuous (see Fig. [T6)) . We note that the first or- 
der phase transition is hardly visible on the caloric curve 
l3{E) whereas it is clearly visible on the magnetization 
curve M{E). A second inhomogeneous phase exists for 

(2) (2) 

E < Ec . It appears precisely at the energy Ec at 
which the homogeneous phase becomes metastable. It 
has a lower entropy S than the homogeneous phase and 
the first inhomogeneous phase and it is unstable. This 
branch is clearly visible on the magnetization curve of 
Fig. 1161 The discussion is similar in the canonical en- 
semble. 

Region 3-c: (i) In MCE, there exists a second order 
phase transition at i?c^' and a first order phase transi- 
tion at Et [i^. (ii) In CE, there exists a second order 

phase transition at /3c^^ and a first order phase transi- 
tion at i3t. For 0.595477 < E < 0.595629, the ensembles 
are inequivalent (see Fig. I15p . However, this concerns a 
strikingly narrow range of energies. 
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FIG. 12: Enlargement of Fig. [TT]in the region of second order 
pliase transition. 
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FIG. 13: Enlargement of Fig. [TT]near the energy Ec where 
a second inhomogeneous phase (unstable) appears. At that 
point, the homogeneous phase becomes metastable. We also 
note, in passing, that the temperatures of the metastable ho- 
mogeneous phase and of the fully stable inhomogeneous phase 
cross each other at some point but this does not signal a 
change of stability. In particular, the homogeneous phase 
remains metastable until the energy Et of first order phase 
transition (see Fig. \T4\i . 



FIG. 15: Enlargement of Fig. [TT] in the region of ensem- 
bles inequivalence. The full line corresponds to the strict 
caloric curve in the MCE and the dashed line to the strict 
caloric curve in the CE. The states between Ei ~ 0.59547 
and E2 — 0.595629 are stable in the microcanonical ensemble 
while they are unstable (i.e. inaccessible) in the canonical en- 
semble. Note that the domain of ensembles inequivalence is 
almost imperceptible. 
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FIG. 16; Magnetization versus energy in Region 3-c showing 
clearly the first and second order phase transitions as well as 
the second (unstable) inhomogeneous branch. 
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FIG. 14: Enlargement of Fig. [TT] in the region of first order 
phase transition. 



E. Region 3-b 

In Figs. [1711^ wc plot the series of equilibria in Region 
3-b corresponding to (/o)i < fo < (/o)2 where (/o)i ~ 
0.109497 and (/o)2 ^ 0.109519 (see Fig. ^. Specifically, 
we consider fo = 0.10950. 

The homogeneous phase exists at any accessible en- 
ergy. It is fully stable for E > unstable for 
e''P < E < E^'^\ metastable for Et < E < eP and 
fully stable ioi E < Et- A first inhomogeneous phase 
exists for E^r^i^ < E < E^^c^ and E < E^ (it does 
not exist between Em^ and i^m'). It is fully stable for 
Et < E < e'^c'^ and metastable for eI^^ < E < Et 
and for E < Ein . Therefore, the microcanonical caloric 
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FIG. 17: Series of equilibria in Region 3-b near the point of 
second order phase transition. 



curve displays a second order phase transition between 
homogeneous and inhomogencous states marked by the 
discontinuity of P'{E) = S"{E) at E = E^c^ and a 
first order phase transition between homogeneous and 
inhomogencous states marked by the discontinuity of 
I3{E) = S'{E) at E = Et. A second inhomogencous 
phase exists for Em^ < E < E^'' and for E < E^ 
(it does not exist between E^"^ and eIu'). It appears 

(2) 

precisely at the energy Ec at which the homogeneous 
phase becomes metastable. This second inhomogencous 
phase is always unstable. The first and second order 
phase transitions are better visible on the magnetization 
curves of Figs. [50] and [HJ The discussion is similar in 
the canonical ensemble. 

Region 3-b: (i) In MCE, there exists a second order 
phase transition at i?c^^ and a first order phase transi- 
tion at Et- (ii) In CE, there exists a second order phase 
transition at /3c^^ and a first order phase transition at /3t. 
As in region 3-c, there exists a tiny region of ensembles 
inequivalence. 
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FIG. 18: Series of equilibria in Region 3-b near the point of 
first order phase transition. We also see the disappearance of 
the inhomogeneous phases for E < . In particular, the 
energy E^ can be interpreted as a spinodal point at which 
the inhomogeneous metastable branch disappears. 
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FIG. 19: Series of equilibria in Region 3-b. We see the re- 
appearance of the inhomogeneous phases for E < Em . The 
energy Ein can also be interpreted as a spinodal point. 



F. Region 3-a 

In Figs. we plot the series of equilibria in Region 

3-a corresponding to (/o)* < /o < (./o)i where (/o)» ~ 
0.10947 and (/o)i ~ 0.109497 (see Fig. [J). Specifically, 
we consider /o = 0.109480. 

The homogeneous phase exists at any accessible en- 
ergy. It is fully stable for E > Ec^\ unstable for 
-B?^ < E < Ei^\ and fully stable for E < E?\ A 
first inhomogeneous phase exists for e'^'^ < E < i^c^-* 
and for E < E^^^ (it does not exist for E^^^ <E < E^^'^). 
It is fully stable for E^^^ < E < Ec^^ and metastable for 

(2) 

E < Em ■ Therefore, the microcanonical caloric curve 
displays two second order phase transitions between ho- 
mogeneous and inhomogencous states marked by the dis- 
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FIG. 20: Magnetization versus energy in Region 3-b. 
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FIG. 21: Enlargement of Fig. [20] in the region of first order 
phase transition. 



FIG. 23: Magnetization versus energy in Region 3-a. It dis- 
plays two microcanonical second order phase transitions at 
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FIG. 22: Series of equilibria in Region 3-a. It displays two 
microcanonical and canonical second order phase transitions 
at Ec^\ Ec'^-' and at /SP' respectively. 
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FIG. 24: Extension of Fig. [23] showing the reappearance of 

(2) 

the inhomogeneous phase for E < Em . 



continuity of I3'{E) = S"{E) at E = E^^'^ and E = Ei'^\ 
The magnetization of the fully stable branch passes from 

M = to M ^ at E = E^^'^ and from M 7^ to 

(2) 

M = at E = Ec , but remains continuous at the 
transition (see Fig. [^5|) . We note that the second order 
phase transitions are hardly visible on the caloric curve 
/?(£') whereas they are clearly visible on the magnetiza- 
tion curve M{E). A second inhomogeneous phase exists 
(2) 

for E < Em and it is unstable. The discussion is similar 
in the canonical ensemble. 

Region 3-a: (i) In MCE, there exists two second order 
phase transitions at iJc^-* and i?c^''; (ii) in CE, there exists 
two second order phase transitions at Pc^^ and l3c^\ The 
ensembles are equivalent. 



0.1075 and (/o)* ~ 0.10947 (see Fig. d]). Specifically, we 
consider /o = 0.10900. 

The homogeneous phase exists at any accessible en- 
ergy and it is fully stable. Therefore, the microcanon- 
ical caloric curve does not display any phase transition 
and is made of homogeneous states. Two inhomogeneous 
phases appear for E < E^, one being metastable and 
the other unstable. The metastable phase has a lower 
entropy than the homogeneous phase and the unstable 
phase has a lower entropy than the metastable phase (see 
Fig. [^5]) . These different phases can also be seen on the 
magnetization (order parameter) curve of Fig. 1271 The 
discussion is similar in the canonical ensemble. 

Region 2: There is no phase transition and the ensem- 
bles are equivalent. 



G. Region 2 H. Region 1 

In Figs. dSjUT] we plot the series of equilibria in Region In Figs. [35] and [221 we plot the series of equilib- 
2 corresponding to (/o)m < /o < (/o)* where (/o)in — ria in Region 1 corresponding to /o < (/o)m where 
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FIG. 25: Series of equilibria in Region 2. There is no phase 
transition but the sudden appearance of a metastable inhomo- 
geneous branch, accompanied by an unstable inhomogeneous 
branch, at i5 = Em- 
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FIG. 26: Entropy versus energy in Region 2. 



(/o)m ~ 0.1075 (see Fig. [T]). Specifically, we consider 
/o = 0.10600. 

The homogeneous phase exists at any accessible energy 
and it is fully stable. There is no inhomogeneous phase. 
Therefore, the microcanonical caloric curve does not dis- 
play any phase transition and is made of homogeneous 
states (see Figs. and YU^ . The discussion is similar in 
the canonical ensemble. 

Region 1: There is no phase transition and the ensem- 
bles are equivalent. 



IV. DISCUSSION 

Let us summarize the different results obtained in the 
previous analysis: 

(i) Decreasing /o, the system successively exhibits one 
second order phase transition (Regions 5 and 4) , one sec- 
ond order and one first order phase transition (Regions 
3-c and 3-b), two second order phase transitions (Region 
3-a), and no phase transition (Regions 2 and 1). 
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FIG. 27: Magnetization versus energy in Region 2. 
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FIG. 28: Series of equihbria in Region 1. 



(ii) There exists a tricritical point corresponding to the 
passage from a first order phase transition to a second 
order phase transition. It is located at ((/o)i, ~ 
(0.109497,0.6059,125). 

(iii) The sudden appearance of two second order 
phase transitions at the turning point ((/o)*, -E^, /?*) ~ 
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FIG. 29: Entropy versus energy curve in Region 1. 
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(0.10947,0.608,118) is sometimes called second order 
azeotropy [27| . 

(iv) For (/o)* < fo < (/o)ci there is a phenomenon 
of p hase reentrance concerning the homogeneous phase 
[rjim ■ As we reduce the energy, the homogeneous phase 
is successively stable, unstable and stable (or metastable) 
again. This phenomenon is basically due to the turning 
point of the energy curve £'c(/o) at /o — {fo)*- It is 
therefore associated with the second order azeotropy. 

(v) For (/o)* < fo < (/o)2, there is a phenomenon of 
phase reentrance concerning the inhomogeneous phase. 
As we reduce the energy, the inhomogeneous phase is 
stable (or metastable), then it disappears, and it finally 
reappears as a metastable state. This phenomenon is 
basically due to the turning point of the energy curve 
E,nifo) at fo = (/o)2- It is located at ((/o)2, -Ea, /32) ^ 
(0.109519,0.603,137). 

(vi) The tricritical point ((/o)i,-E'i) separating first 
and second order phase transitions is located between 
the turning points of the i?c(/o) and i?,n(/o) curves. 

(vii) In Regions 3-c and 3-b, there is a very small zone 
of ensembles inequivalence associated with the first order 
phase transitions. 

In conclusion, the out-of-equilibrium phase transitions 
of the HMF model predicted by the Lynden-Bell theory 
lead to a rich and interesting phase diagram. It is strik- 
ing that everything happens in a very narrow range of 
parameters (/o),„ ~ 0.1075 < fo < {fo)c ^ 0.11253954, 
although fo can take in principle all positive values. A 
similar observation has been made previously in other 
studies of phase transitions in systems with long-range 
interactions [t^, [83 - l86| . The branches corresponding to 
the different phases are very close to each other in the 
series of equilibria j3{E) and in the entropic curves S{E). 
This shows in particular that all the phases have almost 
the same entropy, even the unstable ones. However, the 
branches appear to be well separated in the kinetic caloric 
curve /3kin{E) [82| and in the magnetization curve M{E). 



V. NUMERICAL SIMULATIONS 

The numerical results obtained in [8l| show that the 
phase diagram deduced from the Lynden-Bell theory pre- 
dicts the right phenomenology, even in the small and 
quite complex region located around the tricritical point. 
In this section, we show that, in the region of metasta- 
bility (according to Lynden-Bell's theory), the system 
displays the usual dynamical behavior of systems in a 
metastable state: a "lethargic" evolution during which 
the system is trapped in a given macrostate (metastable) , 
followed by a sudden jump in a different macrostate 
(fully stable). To the best of our knowledge, this is the 
first evidence of metastability during a QSS. Our results 
are obtained from a set of molecular dynamics simula- 
tions carried out after having prepared the system in a 
state belonging to the region of the {fo,E) plane close 
to the first order transition line. For the chosen param- 
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FIG. 30: Magnetization as a function of time for a system of 
N = 20000 particles at energy _E = 0.5901. The four curves 
are different realizations of the same initial distribution with 
fo = 0.1097. The initial condition is a rectangular waterbag 
distribution with magnetization Mo — 0.01413 (see Eq. (28) 
of P). 



eters {fo,E) = (0.1097,0.5901), the homogeneous phase 
is metastable and the inhomogeneous phase is fully sta- 
ble (see Fig. [T|). Following the temporal evolution of 
the magnetization (see Fig. [30)) . one first observes the 
spontaneous relaxation of the system in the unmagne- 
tized phase (metastable), followed by a sudden jump in 
the magnetized state (fully stable). 

Different runs with the same initial distribution show 
that the jump occurs at random times. This is in agree- 
ment with the ordinary behavior of metastable states, 
where the time of the jump depends on the particular 
"realization" . The average time at which the jump oc- 
curs depends on the size of the system and increases 
with A^. This indicates that "collisions" (finite A^ ef- 
fects) play some role in the dynamics. This is relatively 
unexpected since the regime that we are exploring cor- 
responds to the QSS regime where the Vlasov equation 
should be applicable. In all the numerical simulations 
that we have run, the system spontaneously relaxes to- 
wards the metastable state; it never directly reaches the 
fully stable state. The selection of the QSS, among these 
two states, obviously depends on a complicated notion 
of basin of attraction. Our initial condition consists in a 
rectangular waterbag distribution which has a very small 
magnetization Mq ~ 0.01744 since the chosen energy is 
close to the minimum energy state i?mm(/o) — 0.58766 
which is a waterbag distribution with vanishing magne- 
tization [t^, HH . It is likely that this initial condition 
belongs to the basin of attraction of the homogeneous 
metastable state. It is possible that changing the initial 
condition (still with two levels and with the same fo and 
E but no more waterbag) so as to increase the magneti- 
zation Mo may help the system to access directly to the 
inhomogeneous (fully stable) state. On the other hand, 
on the other side of the first order transition line (i.e. 
for smaller values of /o), the homogeneous state becomes 
fully stable while the inhomogeneous state is metastable. 
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In that case, starting from a rectangular waterbag initial 
condition, the system relaxes towards the homogeneous 
state (since it belongs to its basin of attraction) and stays 
there during the whole QSS regime since it is now fully 
stable (numerical simulations not shown) . Additional nu- 
merical simulations are necessary to get a more general 
picture of the QSS metastability for different values of the 
control parameters (/o, E) and different types of initial 
conditions. 



VI. CONCLUSION 

In this paper we have explored the phase diagram ob- 
tained by applying Lynden-Bell's statistical theory to the 
HMF model. We have found a richer phenomenology 
in the {fo,E) plane [zl, HI than in the {Mo,E) plane 
[tsI [79} . We have also explained that the proper exter- 
nal parameter to use in the Lynden-Bell theory is /o, not 
Mf). The choice of the proper control parameters has 
deep consequences on the thermodynamical analysis. 

The HMF model is a system in which the Lynden- 
Bell theory works relatively well (in contrast to astro- 
physical systems for which it was initially devised poj). 
In particular, the phenomenon of phase reentrance that 
is predicted on the basis of this theory [tJ has been 
successfully reproduced in [8l|. Numerical evidence of 
first and second order phase transitions has also been 
given in [8l| in agreement with theory (93] ■ This is 
remarkable because all these interesting features occur 
in a very small region of the phase diagram (typically 
{foU - 0.1075 < fo < (/o)c - 0.11253954). In this 
sense, the Lynden-Bell prediction is not only qualitative 
but in fact extremely accurate! 

There are, however, cases where the Lynden-Bell the- 
ory fails. In the numerical study of [8l|, some discrep- 
ancies with the Lynden-Bell prediction were reported. 
In particular, unmagnetized states are observed in the a 
priori magnetized region leading to a second reentrant 
phase. Inversely, magnetized states are observed in the 
a priori unmagnetized region. On the other hand, the 
Lyndcn-BcU theory cannot explain the region of nega- 
tive kinetic specific heat observed numerically by Antoni 
& Ruffo HH and Latora et al. [tO]- These authors start 
from an initial condition with magnetization Mq = 1 in 
which all the particles are located at 6* = 0. The ini- 
tial distribution function /o is infinite corresponding to 
the dilute (or non degenerate) limit of the Lynden-Bell 
theory in which the predicted QSS coincides with the 
Boltzmann distribution. Now, the results of numerical 
simulations [3^ [70| are inconsistent with the Boltzmann 
(hence Lynden-Bell) distribution in the region of negative 
kinetic specific heats. This means that violent relaxation 
is incomplete [soj and that the system is trapped in a 
stable steady state of the Vlasov equation that is not 
the most mixed (i.e. Lynden-Bell) state [7J|. Recently, 
Chavanis & Campa [t^I have investigated the Vlasov dy- 
namical stability of polytropic (or Tsallis) distributions 



and argued that polytropes with an index close to n = 1 
could provide an explanation of the curious anomalies 
observed in [s^, [l^l- In this work, the polytropic dis- 
tributions are justified by a lack of ergodicity and by 
incomplete relaxation. It would be interesting to extend 
their analysis (restricted so far to initial conditions with 
magnetization Mq = 1) so as to cover a wider range of 
parameters and see whether it can explain similarly the 
anomalies reported in (8l| . 

Very recently, a mathematical "tour de force" has 
been accomplished by Mouhot & Villani [s^] who rig- 
orously proved that systems with long-range interactions 
described by the Vlasov equation possess some asymp- 
totic "stabilization" property in large time, although the 
Vlasov equation is time-reversible. More precisely, they 
show that if a stable steady state of the Vlasov equa- 
tion is slightly perturbed, the perturbation converges in 
a weak sense towards a steady distribution through phase 
mixing without the help of any extra diffusion or ensem- 
ble averaging. This is refered to as nonlinear Landau 
damping. This is a very important work that shades new 
light on the process of phase mixing and, consequently, 
on the nature of QSSs. However, these authors criti- 
cize the Lynden-Bell approach arguing that there is no 
"universal" large time behavior of the solutions of the 
Vlasov equation in terms of just the conservation laws 
and the initial datum. In their words: "This seems to be 
bad news for the statistical theory of the Vlasov equa- 
tion pioneered by Lynden-Bell" . Although it is clear that 
the Lynden-Bell theory has some limitations due to in- 
complete relaxation (lack of mixing/ergodicity) l30| . our 
scries of works related to the HMF model [tJ: l8ll Is^ , 
including the present effort, shows that the Lynden-Bell 
approach is able to make accurate predictions that are 
confirmed by direct numerical simulations. Therefore, 
the Lynden-Bell theory remains a valuable tool even if 
it is difficult to specify its general domain of validity. In 
fact, Mouhot & Villani [s^l do not totally reject this sta- 
tistical approach and point out limitations in the appli- 
cation of their results. In particular, their theory is based 
on smooth functions (which is not the norm in statisti- 
cal theories) and Landau damping is a thin effect which 
might be neglected when it comes to predict the "final" 
state in a "turbulent" situation (which is precisely the 
aim of Lynden-Bell's statistical theory). The subject is 
certainly not closed and should lead again to interesting 
findings and fruitful discussions. 



Appendix A: The ground state 

In this Appendix, we briefly discuss the ground state 
of the Lynden-Bell distribution (analogous to the Fermi- 
Dirac distribution) and its connection with the phase di- 
agram of Fig. [TJ 

For a given value of /o, the minimum energy state cor- 
responds to a Fermi distribution at T = 0, i.e. a (possibly 
spatially inhomogencous) waterbag distribution. Such a 
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FIG. 31: Phase diagram close to the ground state. The up- 
per fuU hues correspond to states that have been actually 
computed. Arrows give the directions towards which these 
branches should tend. The lower line denoted Emin (like 
in Fig. [l| corresponds to the spatially homogeneous wa- 
terbag distribution as explained in [tJ, IsJI- Its energy is 
E-min ~ l/(967r^/Q) ~f 1/2. It is the global energy minimum 
for /o < (/o)i, a local energy minimum for {fo)'t < fo < (/o)c 
and an unstable saddle point for fo > (/o)c- The other lower 
line corresponds to the inhomogeneous waterbag distribution 
as explained in [s^. It starts at fo = (/o)m (corresponding to 
an energy E'^ ~ 0.59473) and tends to S = for fo — >■ +00. 
It is a local energy minimum for (/o)m < fo < (/o)t and a 
global energy minimum for fo > (/o)t- The unstable inhomo- 
geneous waterbag distribution has not been represented (see 
[88 ] for details). 



distribution is equivalent to a polytrope of index n = 1/2 
frit. I ts structure and stability are described in detail 
in [TtI . IssI ]. Here, we only give the final results of the 
analysis (see Figs. [21] and and refer the reader to 
[77I IssI for more details. The global minimum energy 
state (G) is a spatially homogeneous waterbag distribu- 
tion for fo < (/o)t — 0.109579 and a spatially inho- 
mogeneous waterbag distribution for /o > {foYt- On. 
the other hand, the spatially inhomogeneous waterbag 
distribution is a local minimum energy state (L) for 
(/o)J„ ~ 0.1075 < fo < ifoYt and the spatially homo- 
geneous waterbag distribution is a local minimum en- 
ergy state (L) for (foY < fo < (/o)c = 1/(2^^2). For 
fa < (/o)m (spinodal point), no spatially inhomogeneous 
waterbag distribution exists and for fo > (/o)c the spa- 
tially homogeneous waterbag distribution is an unstable 
saddle point of energy. 

We emphasize that the specific form of the initial con- 
dition may constrain the accessible range of energies. 
For example, for a rectangular waterbag initial distribu- 
tion, the minimum accessible energy i?j\//Ar(/o) is strictly 
larger than the ground state Egroundifo) for fo > (/o)t 
(see Fig. [T]). Of course, smaller energies can be achieved 
by other types of initial conditions. 

It is likely that the point (foYt — 0.109579 should coin- 
cide with the transition point {fo)t hi the phase diagram 
of Fig. [T] although we gave a different value (/o)t — 
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FIG. 32: Magnetization of the minimum energy state (ground 
state) as a function of the maximum value of the distribution 
fo (taken from [IHl)- It exhibits a first order phase transition 
at (/o)t- The homogeneous waterbag distribution is the global 
energy minimum for fo < {foYt, a local energy minimum for 
(.fo)t < fo < {fo)c and an unstable saddle point for fo > (/o)c- 
The inhomogeneous waterbag distribution exists only for fo > 
{foYm- It is a local energy minimum for (foYm < fo < {foYt 
(corresponding to M* ~ 0.37 < M < Mt ~ 0.44) and the 
global energy minimum for fo > (/o)t (corresponding to M > 
A4t)- It is an unstable saddle point for (/o)m < fo < (/o)c 
and M < A/.). 



0.10965 in Sec. IIIIl In fact, as we indicated in the cap- 
tion of Fig. [1] the curve £'t(/o) has been continued "by 
hand" for small energies so that the value {fo)t — 0.10965 
is not firmly established and may be incorrect. The 
points that have been actually computed are shown in 
Fig. [211 It is likely that the real curve i?t(/o) tends 
to the point {{foYt, Et{{foYt)) - (0.109579,0.587896). If 
this picture is correct, it implies that the curve Et{fo) is 
multivalued in some range of parameters [(/o)t, {fo)new\ 
(say). Indeed, some of the computed points have val- 
ues of fo larger than {foYt that the curve must turn 
back. This yields an even more complex phase diagram 
with an additional phase reentance. Indeed, decreasing 
the energy in the range [{foYu {fo)new\, the homogeneous 
phase is successively stable, unstable, metastable, stable, 
and metastable again. On the other hand, the inhomo- 
geneous phase is inexistent, stable, metastable and sta- 
ble again. There exists therefore one second order phase 
transition and two first order phase transitions in this 
range of parameters. 

On the other hand, it is likely that the local mini- 
mum energy state for [foYm < fo < {foYt corresponds 
to the minimum accessible energy of the inhomogeneous 
phase while the local minimum energy state for {foYt < 
fo < (./o)c corresponds to the minimum accessible en- 
ergy of the homogeneous phase. As a result, the curve 

Ern {fo) would not connect the point ((/o)m,£^m) — 
(0.1075, 0.59133) as shown in Fig. [T]but rather the point 
{{foYm^E'^/) ~ (0.1075, 0.59473) as shown in Fig. ^ 
Additional numerical simulations would be necessary 
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to ascertain these results but they need to be very ac- proaching the minimum energy, 
curate and we experienced numerical problems when ap- 
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to J as the "free energy" by a slight abuse of language. 

[94] In particular, the thermodynamical specific heat C = 
dE/dT is necessarily positive in CE while it can be pos- 
itive or negative in MCE. 

[95] It can be noted that, in case of incomplete relaxation, 
the initial conditions {E, Mq^') and {E, M^^) can lead to 
different results, although the Lynden-Bell theory leads 
to the same prediction. This could be interesting to check. 

[96] A similar situation has been observed in [tJ for poly- 
tropic distributions. 

[97] We should be careful, however, that Staniscia et al. 
[sH do not determine the detailed distribution functions 
f{6, v) of the inhomogeneous states. It is possible that the 
Lynden-Bell theory predicts the correct phase transitions 
in terms of the magnetization although the distribution 
functions are not exactly given by Eq. ([6J. This would be 
interesting to check in future works. 



